arXiv: 1505.02787v2 [astro-ph.CO] 9 Sep 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 10 September 2015 (MN DTeX style file v2.2) 


The clustering of Lyman alpha emitters at z « 7: 
implications for reionization and host halo masses 

Emanuele Sobacchi 1 * & Andrei Mesinger 1 

1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy 


10 September 2015 


ABSTRACT 

The Lya line of high-redshift galaxies has emerged as a powerful probe of both early galaxy 
evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-held sur¬ 
vey with the Subaru Hyper Suprime-Cam (HSC), we study the angular correlation function 
(ACF) of narrow-band selected, z ~ 7 Lyct emitting galaxies (LAEs). The clustering of LAEs 
is determined by both: (i) their typical host halo masses, M h ; and (ii) the absorption due to 
a patchy EoR, characterized by an average neutral fraction of the intergalactic medium, xhi- 
We bracket the allowed LAE ACFs by exploring extreme scenarios for both the intrinsic Lya 
emission and the large-scale pattern (i.e. morphology) of cosmic ionized patches in physical 
EoR models. Independent of the EoR morphology, current z ~ 7 ACF measurements con¬ 
strain x m < 0.5 (l-cr). We also find that the low values of the currently-observed ACF imply 
that LAEs are hosted by relatively small dark matter halos: M h < 10 1 o Mq, with correspond¬ 
ing duty cycles of < few per cent. These values are over an order of magnitude lower than 
the analogous ones for color-selected, Lyman break galaxies, suggesting that z ~ 7 narrow- 
band LAEs searches are preferentially selecting young, star-burst galaxies, residing in less 
massive halos. The upcoming Ultra Deep campaign with the HSC will significantly improve 
constraints on both the EoR and LAE host halos. 

Key words: cosmology: theory - early Universe - dark ages, reionization, first stars - galax¬ 
ies: formation - high-redshift - evolution 


1 INTRODUCTION 

Although it is the last major phase change in the history of our Uni¬ 
verse, the epoch of reionization (EoR) remains poorly explored. 
The EoR is expected to be fairly extended and patchy, with the 
growth of cosmic HII regions sourced by the birth of an increasing 
number of early galaxies. Understanding this complex process re¬ 
quires accurate statistics. Luckily, efforts are underway to increase 
the sample of high-redshift (z > 6) objects serving as EoR probes, 
such as quasars and galaxies, hopefully resulting in an unambigu¬ 
ous EoR detection. 

Some of the most important of these efforts are focused on Ly¬ 
man alpha emitters (LAEs). Due to the significant Lya damping 
wing opacity of the neutral intergalactic medium (IGM), Lyman 
alpha emission from galaxies located close to the edges of cos¬ 
mic HII regions will be strongly attenuated. Conversely, since their 
Lya photons can redshift even further away from resonance before 
reaching the neutral IGM, galaxies far from these edges will be less 
affected and therefore easier to observe. Thus the EoR will impact 
both the observed number of LAEs as well as their clustering prop¬ 
erties. 

* email: emanuele.sobacchi@sns.it 


The main difficulty in using LAEs as an EoR diagnostic lies in 
the fact that we do not understand their intrinsic properties, such as 
their Lya luminosity, the profile of the Lya lines, and the absorp¬ 
tion within the local circumgalactic medium (CGM). These intrin¬ 
sic properties can be degenerate with an EoR signature. Lacking 
this a priori information on intrinsic properties, most modern stud¬ 
ies focus on redshift evolution, with the implicit assumption that 
most of the evolution is due to the EoR: a change in the volume- 
average neutral fraction of the IGM, aim. Indeed, both the fraction 
of drop-out galaxies with strong Lya emission and narrow-band 
LAE surveys show a drop in these populations at 2 > 6, counter- 


ing empirical trends from lower redshifts (e.g. Ouchi et al. 

2010; 

Fontana et al. 2010 Stark et al. 2010; Kashikawa et al. 2011; 

Caru- 

ana et al. 2012, Treu et al. 2013, Konno et al. 2014; Schenker et al. 

2014; Caruana et al. 2014; Pentericci et al. 2014, Cassata 

et al. 

2015). However, small-number statistics, as well as ignorance of 


the evolution in intrinsic Lya emission and its relation to contin¬ 
uum selected galaxies prevents a robust claim of an EoR detection 


(e.g. Dijkstra & Wyithe|2010 

Dayal et al.|2011 Bolton & Haehneltj 

2013ljensen et al. 2013 Dij 

cstra et al. 2014 Taylor & Lidz 2014 

Mesinger et al. 2015). 


An alternative EoR diagnostic is the observed clustering of 
LAEs. Clustering could be a more powerful EoR probe than pure 
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number counts, for several reasons. Firstly, the intrinsic (not modi¬ 
fied by the IGM) clustering of LAEs is determined almost entirely 
by their host halos. The correlation functions of dark matter (DM) 
halos are well understood and only evolve by factors of ~few over 
the range of possible host halo masses; this is in contrast to un¬ 
certainties in baryonic physics (e.g. the Lya escape fraction, star- 
formation rates) which can have an order unity impact on the ob¬ 
served number density. Secondly, the galactic environment (includ¬ 
ing self-shielded systems and accretion flows), while possibly hav¬ 
ing a large (degenerate with that of the EoR) imprint in the number 
density and its evolution (e.g. |Dijkstra & Wyithe|2010||Finkelstein| 
|et al.|20T2{ |Bolton & Haehnelt|2013[|Mesinger et al ,|20 1 5f , will 

have a much weaker impact on the the observed clustering signa¬ 
ture on the large (~tens Mpc) scales relev ant for the EoR (e.g.|Cro-| 
|ciani et al.|201 l|[Wyithe & Dijkstra|201 l||Mesinger et al.|2015|>. 

The above means that the uncertainties associated with the in¬ 
trinsic clustering signature of LAEs are much smaller than those of 
their intrinsic number densities. This could allow us to constrain the 
EoR properties from a just a single, high -2 clustering measurement 
(e.g. Furlanetto et al.|2006[[M cQuinn et al.|2007fr. Indeed, using a 
counts-in-cell statistic, [Mesinger & Furlanetto[ l |2008bj > showed that 
the LAE clustering induced by patchy reionization is so strong dur¬ 
ing the first half of the EoR, that it exceeds even the worst case sce¬ 
nario of LAEs being hosted by the most massive, most intrinsically- 
clustered halos. In other words, the EoR-induced clustering would 
have an unambiguous imprint during the first stages of reioniza¬ 
tion, which cannot be replicated by even ad-hoc models for relat¬ 
ing intrinsic Lya luminosities to their host halos L“ tr Mu. 
Such a measurement would provide the much-needed “smoking- 
gun” proof for reionization. 

Unfortunately, we are unlikely to have a measurement of the Lya 
clustering at the redshifts high enough to probe the first half of 
reionization (which best-fit [Planck Collaboration et al.|2015[ mod- 
els place at z > 9, albeit with large error bars). Indeed, preliminary 
estimates suggest that the Ultra Deep campaign with the HSC will 
only capture tens of LAEs at z ~ 7.3 (compared to > 1000 at 
2 « 6.6; M. Ouchi, private communication). Hence, constraints 
using LAE clustering during the second half of reionization (as ex¬ 
pected to be the case at z « 6.6) must explore the degeneracy 
between which halos host LAEs (i.e. intrinsic clustering) and the 
patchy reionization-induced clustering. 

In this study we compute the angular correlation functions 
(ACFs) of LAEs at z ~ 7, using several EoR models as well as 
prescriptions for assigning intrinsic Lya luminosities to host halos 
(L^ tT o Mu). Our work is similar to the recent predictions by 
Jensen et al.~| (2014| ; however, given the large modeling uncertain¬ 
ties at high redshifts, here we adopt a more systematic approach 
by using extreme scenarios both for EoR morphologies as well as 
the L™ r -O' Mu relation in order to bracket the allowed parameter 
space. Moreover, our framework uses the latest EoR simulations 
from |Sobacchi & Mesinger] ( |2014) >, which include state-of-the-art 
sub-grid models of important physical processes like inhomoge¬ 
neous recombinations inside small-scale structures. As this work 
was nearing completion, a similar effort was also made by [HutterJ 
[et al.[(2015| ). Our work also differs from this analysis, which de¬ 
rives EoR constraints using one specific LAE model and one spe¬ 
cific reionization scenario driven by z = 6.6 galaxies. 

This study is motivated by the upcoming 2 = 6.6 Subaru survey 
with the Hyper Suprime-Cam (HSC). Current Subaru clustering 
measurements are suggestive of a mostly ionized universe at 2 « 7 
( [McQuinn et al.|2007HQuchi et al.|2010[|Jensen et al.|20 1 3} . The 

upcoming HSC survey will result in up to a dramatic 30-fold in¬ 


crease in survey volume (~30 deg 2 for the deep survey; M. Ouchi, 
private communication), finally allowing us to statistically sample 
the EoR morphology (Taylor & Lidz|2014|| Jensen et al. |20 1 4] >. 

This paper is organized as follows: in Section [2] we present our 
models for reionization and LAEs; in Section [3] we see how to re¬ 
cover robust limits on xhi from clustering measurements; in Sec¬ 
tion [4] we present our conclusions. Throughout we assume a flat 
ACDM cosmology with parameters (D m , Da, Db, h, og, n) = 
(0.28, 0.72, 0.046, 0.70, 0.82, 0.96), as measured by the Wilkinson 
Microwave Anisotropy Probe (WMAP; [Hinshaw et al.[20 13), c on- 
sistent with recent results from the Planck satellite (Planck Collab-| 
|oratio n et al.|2015} . Unless stated otherwise, we quote all quantities 
in comoving units. 


2 MODEL 

The observed Lya luminosity of a given galaxy, residing in a DM 
halo of total mass, Mu, can be written as: 

L a = LL ntr e- TIGM > L™ in , (1) 

where LJj 1 T (Mh) is the intrinsic (before passing through the IGM) 
Lya luminosity, tigm(xhi) is the IGM optical depth along a given 
sightline, and L™ ln (z) is the sensitivity threshold of the observa¬ 
tion (e.g. = 2.5 x 10 42 erg s _1 for the upcoming HSC 

Ultra Deep survey at 2 = 6.6). Thus, modeling the clustering of 
observed LAEs requires two components: (i) a scheme for assign¬ 
ing intrinsic Lya emission to host DM halos, and (ii) 

large-scale reionization simulations for determining the IGM opac¬ 
ity during the EoR, tigm^hi). Given that both of these are poorly- 
constrained, we adopt a systematic approach and explore extreme 
models for both (i) and (ii). Below we describe each in turn. 


2.1 The intrinsic emission of LAEs 


The intrinsic Lya emission of high -2 galaxies depends on compli¬ 
cated physical processes inside the interstellar medium (ISM) and 
circumgalactic medium (CGM), and cannot be computed from first 
principles. Instead, intrinsic luminosities can be predicted empiri¬ 
cally using observations of the Lya equivalent widths at lower red¬ 
shifts (where the LAE samples are sufficiently large), which relate 
the Lya luminosity to the UV continuum magnitude through some 
probability distribution function (see, e.g. jDijkstra & Wyithe|2012| 
and references therein). The UV continuum magnitude can then be 
related to the DM halo mass using standard abundance matching, 
and extrapolating down the faint end slope, under the_assumption 


Kuhlen & Faucher- 


of a mass-independent duty cycle, fdut 
Giguere 20121. This procedure however is highly uncertain, rely 
ing on empirical, poorly-constrained extrapolations in both redshift 
and magnitude, as well as implicit assumptions about the mass and 
redshift scalings of distributions and duty cycles. 

Here we adopt a different approach. We explore three paramet¬ 
ric prescriptions relating the Lya luminosities to halo masses. Al¬ 
though these models match current observational constraints, they 
are not intended to be overly realistic. Rather they are intended to 
bracket the allowed range of the clustering signal. Consistent with 


1 Here we define a duty cycle as the fraction of halos of a given mass which 
host detectable LAEs, i.e. with L a > 2.5 X 10 42 erg s —1 for our fiducial 
mock survey. 
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z ~ 4 observations (e.g. |Gronke et al.|2Q15^ , we assume the Lya 
luminosity increases with halo mass: 


L intr _ 7- n 

rv — J- / n 


M h 

M™ ir 


( 2 ) 


where \ is a random variable (x = 1 with probability /duty and 
X = 0 otherwise), which accounts for the expected bursty star for¬ 
mation inside high-z DM halos (e.g. Lee et al. 2009, Wyithe et al. 
|2014| >. The normalization of the intrinsic luminosity in eq. |2]( is 
governed by the halo mass corresponding to the detection thresh¬ 
old, M™ ln , and the shape of the LF can be parametrized with the 
power law index, /3. Our fiducial models assume 6 = 1, though we 
also compute correlation functions with (3 = 2/3 without notable 
changes to our conclusions. 

Our systematic approach therefore makes use of three “tuning 
knobs” for the intrinsic luminosity of LAEs: (i) the Lya duty cycle, 
/duty, which shifts the LAE LFs up and down; (ii) the normal¬ 
ization, M™ ln , which shifts the LAE LFs left and right, and (iii) 
the power-law scaling with mass, /3, which flattens or steepens the 
LAE LFs. For a given value of (ii) and (iii), changing (i) does not 
impact the LAE clustering properties, only the observed number 
density of LAEs. In our analysis of the z = 6.6 clustering, we 
vary /duty with the neutral fraction of the IGM during the EoR, 
in order to fix the observed number density of LAEs to the already- 
constrained value of iilae( z = 6.6) = 4 . ll 9 ' 9 x 10~ 4 Mpc" 3 (for 


2.5 x 10 42 ergs" 1 : 


~ - LA-n a ^ ^ 

Ouchi et al. 2010). Therefore, 


a given 


value of M""' : requires a higher value of /duty earlier in reioniza¬ 
tion, in order to compensate for the added attenuation of the IGM^ 
Instead, (ii) and (iii) do impact the intrinsic LAE clustering. 
Their combined impact can be better parametrized in term of the 
luminosity-weighted average host halo mass: 


Mh = 


M?n(M h )dM h 

f+Zn n(Mh)dMh 


(3) 


where n(AIh) is the halo mass function (HMF) at z = 6.6 (e.g. 
|Jenkins et al.|200l] ). The HMF is quite steep over the relevant range, 
so that Mh is at most a factor of few higher than M™ ln for rea¬ 
sonable values of ft ~ 1, with the faintest, most abundant LAEs 
dominating the intrinsic clustering. Below, we keep /3 — 1 and 
vary the normalization in order to change the intrinsic clustering 
properties, as parametrized by the luminosity-weighted mean halo 
mass. Although they do not impact our clustering conclusions, we 
do show some models with different values of /3, illustrating that 
more accurate LFs could constrain this scaling. 

We take three models for Mh, spanning the allowed range: 

• Most massive halos, M h ~2x 10 11 Mg: this relation max¬ 
imizes the intrinsic clustering of LAEs by assuming that the ob¬ 
served LAEs are hosted in the most massive, most intrinsically 
clustered DM halos. This is achieved by using a duty cycle of 
unity, /duty = 1, which results in an average halo mass of Mh s=s 
1.8 x 1O 11 M 0 , when normalized to the observed LAE number 
density. 

• Least massive halos, Mh ~3x 10 9 Mq\ this relation is cho¬ 
sen to minimize the intrinsic clustering of LAEs. LAE populate 


2 The exception to this procedure is the “most massive halo” model below, 
which assumes a duty cycle of unity. As the duty cycle cannot be further 
increased, for this model we allow lower observed number densities during 
the EoR. However, as we shall see below, this extreme model is already 
ruled out. 



Figure 1. LAE luminosity functions at redshift z = 6.6; dots/squares rep¬ 
resent observations by |Ouchi et al.| (2010l/Matth ee et a i~~] (21) 1 5) . With the 
solid, dashed and dotted lines we show the simulated luminosity functions at 
Shi = 0 with different choices of the average mass of halos hosting LAEs 
(Mh = 2 X 10 11 , 2 X 10 10 and 3 X 10 9 Mq respectively, and = 1), 
which show a good agreement with the data. For comparison we show the 
LAE luminosity function assuming Mh = 3 X 10 9 Mq and ft = 2/3. 


DM halos down to masses of M™ ln = 10 9 Mq. This is an extreme 
value, corresponding to the conversion of all halo baryons into stars 
inside a single dynamical timescale, with all Lya photons escaping 
the galaxyFj Matching the observed LAE number densities in this 
model requires a duty cycle /duty ~ 0.001, assuming a mostly- 
ionized Universe. 

• Moderately massive halos, Mh ~ 2 x 10 10 Mq: this rela¬ 
tion is chosen as an intermediate between the two extremes above, 
and is in approximate agreement with the previous best-fit values 
obtained from the clustering of LAEs at z ~6-7 (e.g. |Ouchi et al.| 
|2010j ). The corresponding duty cycle is /duty ~ 0.02, assuming a 
mostly-ionized Universe. 


1 The minimum mass of the halos hosting visible LAEs (with L/ lt,r 
~ 2.5 X 10 42 ergs -1 can be crudely estimated by the following ar¬ 
gument (c.f. [Mesinger & burlanctto|2008b| : assuming that about 2/3 of 
the ionizing photons absorbed within the galaxy are converted into Lya 
photons ( Osterbrock| [l989|, one can write the conversion as I' r = 
0.67 X hv a (1 — /esc) A e -,7V res , where v a is the rest-frame Lya fre¬ 
quency, fesc is the escape fraction of ionizing photons, p * is the star 
formation rate (SFR), e 7 is the ionizing photon efficiency per stellar 
mass and T 7 , res is the fraction of Lya photons which escape from the 
galaxy without getting resonantly absorbed. Assuming that galaxies con¬ 
vert a fraction, /*, of their gas into stars over some mean time-scale, 
t„, and that f esc -C 1, one can write the above relation as: L/ lt,r ~ 
2xl0~ 12 ( ) Mh erg s“ 1 . We obtain our most extreme model 

for LJJ ltr (Mn) by maximizing £ 7 /„T 7 , res /t*: £ 7 = 6 X 10 6 ° ionizing 
photons Mfi 1 ( corresponding to a PopII IMF with Z = 0.04 X from 
|Schaerer|2002) , /*T 7 , res = 1, t* ~ 200 Myr corresponding to the dy¬ 
namical time at z ~ 7. These extreme assumptions result in the relation 
L)/ tr ~ 2 X 10 42 (Mu/10 9 Mq) erg s" 1 , corresponding to a minimum 
mass-scale ~ 10 9 Mq for the faintest LAEs in the Subaru Ultra Deep sur¬ 
veys. Nevertheless, this estimate is very rough, and it is theoretically possi¬ 
ble to detect extremely young star-bursts in any halo massive enough to host 
a galaxy (above the atomic cooling threshold). We explore this even-more- 
extreme possibility in the Appendix (Fig. |A3) , showing that our conclusions 
remain unchanged, since ACFs are very weak functions of halo mass this 
far below the knee of the HMFs. 
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Figure 2. Angular correlation func tions of L AEs at z= 6 . 6 ; black and grey 
dots represent observations by Ouchi et al.j|2010l for LAEs with L“ lm ~ 
4.0 X 10 42 erg s _1 and L™ ln ~ 2.5 X 10 42 erg s _1 respectively. With 
the solid, dashed, dotted lines we show the simulated correlation functions 
at Shi = 0 when Mh = 2 X 10 11 , 2 X 10 10 and 3 X 10® Mg respectively. 


2.1.1 Instrinsic line profiles 


The intrinsic profile of the Lya line is set by radiative transfer 
trough the ISM and CGM, and depends on dust, DLAs, geometry 
and gas kinematics (e.g. see the recent review by |Di jkstra|2014| . 
All of these processes are highly uncertain in high -2 galaxies. Here 
we merely evaluate the IGM absorption at a fixed velocity off¬ 
set from the systemic r edshift, chosen in our fiducial model to be 
A v =200 km s -1 (e.g. Shibuya et al.|2014 Stark et al.|2014 


So- 


bral et al.|2015 i. We also show results using Av =100 km s in 


Fig. Al, which serve to highlight that, unlike EoR constraints based 
on the (evolution of the) LAE number density, our large-scale clus¬ 
tering constraints are extremely insensitive to the chosen line pro¬ 
file. This is because we do not discriminate EoR models based on 
the observed LAE number density, tilae(2 = 6.6), as this quantity 
depends on both intrinsic and EoR properties. As discussed above, 
we instead evaluate clustering at a fixed number density (see also 
Jensen et al.|2014) , effectively adjusting our free parameter /duty 
to compensate for the stronger/weaker IGM absorption due to a 
smaller/larger systemic velocity offset]^] 


2.1.2 Comparison to current z = 6.6 observations, assuming an 
ionized Universe 

We construct our mock survey by slicing the 2 = 6.6 halo field 
of our simulations (see §2.2| l into slabs of width ~ 40 Mpc, cor¬ 
responding to the redshift uncertainty, Az ~ 0.1, for the narrow- 
band LAE surveys fOuchi et al.|201(j) . Each halo is assigned an 
intrinsic Lya luminosity according to eq. Q- and the IGM opacity 
is computed by stepping through the simulation from the halo cen¬ 
ter out to a distance of ~ 200 Mpc along the chosen line of sight 


Figure 3. Evolution of the volume-weighted global neutral fraction, xm 
for the SMALL.HII (solid) and LARGE.HII (dashed) models; the 
LARGE-HII model corresponds to EoR galaxies hosted only by the most 
biased halos, whereas the SMALL_HII EoR model is driven by less biased, 
smaller-mass galaxies at the atomic cooling threshold. This figure serves to 
illustrate physically-motivated values of 5 ri and its late-time evolution. In 
keeping with the systematic approach of this study however, we use the 
neutral fraction fields from various values of 5 hi when computing LAE 
clustering at 2: « 7 (i.e. we do not assume a prior on Sri (z ~ 7)). 
We note that the SMALLJHLII (solid) and LARGEJHII runs have r e = 
0.07 and 0.05 respectively, both within the recent l-cr constraints published 
by Planck Collaboration 2015. We also note on the plot the strict, model- 
independent upper limits around z ~ 6 from the dark fraction in QSO 
spectra (McGreer et al. 2015). 


/■ 


cdt 


tigm = / dz-^nnctaXm 


(4) 


where tih is the local hydrogen number density and a a is the 
Lya cross-section (taking into account proper motion of the cos¬ 
mic gas). 

Before discussing our reionization simulations, we first show 
LAE properties assuming xm = 0 at 2 = 6.6^] All the quanti¬ 
ties quoted below are calculated taking the average between all of 
the mock-survey slabs. 

In Figure|T]we show the luminosity function of the LAEs at red¬ 
shift 2 = 6.6. Dots represent observations by |Ouchi et al. (|2010|>, 
while diamonds represent the recent determination by Matthee| 
jet al.|{2015| ). The curves correspond to our simulated luminosity 
functions. With the lines we show our simulated luminosity func¬ 
tions assuming xm = 0. The assumption that L™ tT oc M^Ci.e. 
f) = 1) is at odds with the recent observations byJMatthee et al.| 
< |2015) > at the bright end, i/LAEs are hosted by the most massive 
halos (/duty ~ 1)- Although one could flatten the LFs by assuming 
a weaker dependence on halo mass (lower /3), we show below that 
this extreme most massive halos model is ruled out with clustering 
measurements. Although here we assume a post-reionization Uni¬ 
verse, we note that reionization should not appreciably change the 


4 This approach is different from the recent work of Hutter et al. f 2015 1 , 
who instead fix /duty = L and then use a line-centered Gaussian Lya 
profile in order to constrain how ttlae changes during the EoR (which is 
in their model driven by the 2 = 6.6 galaxies themselves). With these added 
assumptions, they were able to obtain much tighter limits on Sri than we 
quote below. 


5 As discussed below, our EoR simulations have a sub-grid model for the 
percent-level residual fraction of HI inside ionized cells. This means that 
tIGM i s not exactly zero post-reionization, and is slightly different in our 
two EoR models, owing to their different mass-weighted neutral fractions 
post-reionization (~ 1% and ~ 3%). These minor differences do not im¬ 
pact our conclusions, as any significant level of absorption post-reionization 
by proximate DLAs (e.g. jMesinger et al.|2Q15) , is by definition taken into 
account (on average) by our L)f ltr «->■ relation. 
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shape of the observed LFs, only its amplitude (e.g. Furlanetto et al. 

2006 

McQuinn et al. 2007, Mesinger & Furlanetto 2008b, Dayal 

et al. 

2011 ). 


In Figure[2]we show the LAE ACFs at z = 6 . 6 , again assum- 
ing aim = 0. Black and grey dots represent observations by |Ouchi1 
et al.||20 lof for LAEs with Z/” ln — 4.0 x 10 42 erg s _1 and L™ m ~ 
2.5 x 10 32 erg s _1 respectively. With the solid, dashed and dotted 
lines we show the simulated ACF for our three intrinsic luminosity 
relations: Mh = 2 x 10 11 , 2 x 10 10 and 3 x 1O 9 M 0 respectively. 
We point out that, since LAEs with L ™ tr > 4.0 x 10 42 erg s _1 
provide only upper limits on the ACF, they can not be used to dis¬ 
criminate models with different Mh; instead LAEs with a lower lu¬ 
minosity L™ tT > 2.5 x 10 42 erg s _1 provide useful constraints on 
Mh, and exclude the most extreme Mh = 2 x 10 11 Mq, as already 
noted by |Ouchi et al.|| 20 To| ). Here we are assuming aim = 0 ; how¬ 
ever, a higher neutral fraction increases the ACF at the same Mh, 
strengthening this conclusion. 


2.2 Large-scale simulations of reionization 


To model cosmological reionization, we use 21CMFASTv2 


(Sobacchi & Mesingerj 2013al 


, a modified version of 

the publicly-available code 21CMFAS1 

! 

(Mesinger & Furlanetto 


of the code below; interested readers are encouraged to read the 
introductory papers. 

We generate the IGM density and source fields by: (i) creating 
a 3D Monte Carlo realization of the linear density field in a box 
with sides L = 500 Mpc and N = 2500 4 grid cells; (ii) evolving 
the density field using the Zeldovich approximation (|Zel’Dovich| 


1970), and smoothing onto a lower-resolution N = 50m"grid; (iii) 
using excursion-set theory i Press & Schechter 1974; Lacey & Cole| 


|1993| |Sheth & Tormen| 1999 [ on the evolved density field to com¬ 

pute the fraction of matter collapsed in halos bigger than a threshold 
halo mass, assumed to be large enough to host star-forming galaxies 
(see below). The Lagrangian positions of z « 7 halos are identi¬ 
fied according to the standard spherical collapse model, and subse¬ 
quently displaced using the Zeldovich approximatiorQ (Mesinger 
|& Furlanetto|200~7| . 

The ionization field is computed by comparing the cumulative 
number of ionizing photons to the number of baryons plus cumu¬ 
lative recombinations, in spherical regions of decreasing radius R 
(i.e. following the excursion-set approach of |Furlanetto et al.|2004) . 
Specifically, a cell located at spatial position and redshift, (x, z), is 
flagged as ionized if: 


^/coll(x, Z, R, AT m in) 2 1 -f- ttrec (x, Z , 7?) 
where the ionizing efficiency can be written as: 


« = 30 <3Mo)(IS 


/* 

0.05 


(5) 


( 6 ) 


and N 1 is the number of ionizing photons per stellar baryon, f esc 
is the fraction of UV ionizing photons that escape into the IGM, /« 
is the fraction of galactic gas in stars, /b is the fraction of baryons 
inside the galaxy with respect to the cosmic mean S 2 b/H m , and 
/coil (x, z, R, M m i n ) is the fraction of collapsed matter inside a 


6 http://homepage.sns.it/mesinger/Sim 

' We confirm that using second order Lagrangian perturbation theory only 
changes the resulting za7 correlation functions by < 2-5%, on the scales 
considered in this work. 


sphere of radius R residing in halos larger than M m i n -|^] Starting 
from the box size, the smoothing scale is decreased (averaging also 
over M m i n (x,z) and n re c(x,z)), and the criterion in eq. |5j is re¬ 
evaluated to determine if a cell is ionized or not. For the remaining 
neutral cells, the partial ionizations from sub-grid sources are ac¬ 
counted for by setting the ionized fraction of the remaining neutral 


Cells tO £/coll(x, Z , RceWi -^Tmin)/(1 H - ^rec) ( 

Zahn et al. 2011 1 . 

As detailed in Sobacchi & Mesinger] (2014 

i, the main improve 

ment of 21 CMFASTv2 is for each cell to keep track of its local val 


ues of M m i n (x,z) and n rec (x,z), computed according to the cell's 
density and ionization history. These values depend on density and 
ionization structure on scales much smaller than the cell sizes of 
EoR simulations, and therefore must be accounted for in a sub-grid 
fashion. 

The minimum halo mass hosting star-forming galaxies can be 
expressed as: 


-Mmin — max [M coo l, .Merit] 5 


(7) 


where M coo \{z) corresponds to the atomic cooling threshold of ha¬ 
los with virial temperatures > 10 4 K, and M cr it(J 2 i, zin, z) cor¬ 
responds to a feedback scale below which star-formation is very 
inefficient in regions exposed to a UV background with intensity 
J 21 since zin (see below). Likewise, the local cumulative number 
of recombinations per baryon is computed according to: 

/ \ f dtlrec dt , 

n rec (x,z)= -ar Tz dz, W 

J z IN 

with 


dn It 


dt 


r+oo 

= Jo 


Py (A, z) AtihOb [1 — xm (A)]" dA . 


(9) 

where Py (A ,z) is the full (non-linear, sub-grid) distribution of 
overdensities: A = nn/hn (e.g. 


Miralda-Escude et al. 


2000 


Sobacchi & Mesinger 2014), and xhi (A) is equilibrium neutral 


fraction computed using the empirical self-shielding prescription 
of |Rahmati et al.| ( |2013| . 


2.2.1 EoR Runs 


Given the current uncertainties on the reionization history and mor¬ 
phology, we compare the results of two different runs, which should 
span the range of possible EoR morphologies: 


• SMALL_HII: corresponding to the fiducial simulation 
in Sobacchi & Mesinger l |2014| , this run ignores inter¬ 
nal feedbacij^] with only photo-heating feedback suppress¬ 
ing star formation. Specifically, we take £ = 30 and 

M crit = M 0 J 2 i“ (At £) 6 [l - ( ' 1 ^ 7 ) 1 , with best-fit val 


ues from Sobacchi & Me singer| < |2013b) : (Mo,a,b,c,d) = 
(2.8x10 Mq, 0.17,—2.1,2.0,2.5). The negative impact of 
photo-heating feedback on the baryon budget at z > 6 is mod¬ 
est on its own (e.g. |Okamoto et al.||2008| |Sobacchi & Mesinger| 


8 Note that the minimum halo mass of cosmic reionization sources, M m [ n , 
is a different quantity than the halo mass corresponding to the faintest z = 
6.6 LAEs discussed above, M™ in . 

9 Supemovae can blow out gas from small potential wells, thereby dra¬ 
matically reducing subsequent star formation inside the smallest halos (e.g. 
|Springel & Hernquist|2003[|Pawlik & Schaye|2009| |Finlator et al.|2011fr . 

The efficiency of such internal feedback and what precisely is the associ¬ 
ated suppression scale, depend on the details of SNe explosions, cooling, 
and ISM structure, and are extremely uncertain at high-redshifts. 
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Figure 4. Power spectra of the ionization fields at x h e = 0.75, 0.50, 0.25 (left to right panels) for the SMALL_HII (solid) and LARGE.HII (dashed) EoR 
models. 


|2013a|[Noh & McQuinn||20T4l |Pawlik et al.||2015| >. However, its 
importance is dramatically increased when inhomogeneous recom¬ 
binations are taken into account ( jSobacchi & Mesinger}2014) , as 
both processes preferentially slow the expansion of the largest HII 
regions, which started growing early-on, thus allowing sufficient 
time for both feedback processes to take their toll. The resulting 
EoR models are extended and characterized by small cosmic HII 
patches. 

• LARGE-HII: this run assumes that star-formation is efficient 
only in the rarest, most-massive, most-biased halos. Such a model 
could be motivated by efficient internal feedback inside small mass 
galaxies (e.g. Springel & Hernquist 2003, Raicevic et al.||201 I 


|Kim et aTT||20 1 3} . We adopt the extreme values of £ = 50 and 

M cr it (T v ir = 10 5 K). Such a high threshold virial temperature ap¬ 
proximately corresponds to (a factor of ~2 lower than) that inferred 
from the faintest observed high-z galaxies, as computed from abun¬ 
dance matching (e.g. |Kuhlen & Faucher-Giguere||2012| |Greig &| 
|Mesinger|20 1 5| . Driven by rapidly-forming galaxies on the high- 
mass tail of the mass function, the resulting EoR models are rapid 
and characterized by large cosmic HII patches. 

In Figure [3] we compare the evolution of the global volume- 
averaged neutral fraction aim in these two different reionization 
models. This serves only to illustrate the plausible range of EoR 
evolutions: having a mostly neutral Universe at 2 = 6.6 is diffi¬ 
cult to achieve with physical EoR models (see also, e.g. |Lidz et al.| 
|2007| . However, in keeping with the systematic nature of our study, 
we do not restrict the allowed values of a im at z ~ 7, instead com¬ 
puting the resulting LAE clustering with ionization maps at any 
*hi(z) superimposing them on the z ~ 7 halo field (e.g. McQuinn 
|et al.|2007| [Mesinger & Furlanetto|2008b| >. 

Therefore, much more relevant to this work is the EoR morphol¬ 
ogy of the two models, which is quantified in Fig. [4] We plot the 
spherically-averaged power spectra, A)^. = k 3 / (2-iHV) (|<5 xx | 2 )fc, 
with <5 XX = xni/xm — 1 at xm = 0.75, 0.5, 0.25. As expected, the 
resulting EoR morphology is very different in the two models, with 
LARGE_HII having ~2-5 times more power than SMALL_HII 
on large (k <0.1 Mpc” 1 ) scales. 

Finally, we note that, unlike recent work on the Lya frac¬ 
tion in Mesinger et al.| ( |2015^ , we do not resolve the opacity of 
self-shielded systems. Although originally thought to be impor¬ 
tant in governing the Lya transmission at high-redshifts (Bolton &| 
| Haehnelt|20 1 3| , their impact is dramatically diminished with more 
realistic self-shielding prescriptions (Mesinger et al. 2015 see also 
|Keating et aL| 20f4i. Moreover, due to their steep damping wing 


have a strong imprint on the Lya line. This means that their poten¬ 
tial impact on the detectability of LAEs can be mostly subsumed 
within the uncertainties in modeling the L™ tT -O' Mh relation dis¬ 
cussed above, and is highly unlikely to have a strong imprint on the 
clustering signature on the ~10 Mpc scales relevant for the EoR 
(e.g. Crociani et al.}20lT| . 


3 RESULTS 

3.1 Mock LAE maps 

In Figure [5] we show the distribution of the observed LAEs (black 
points) in our mock surveys. Columns correspond to increasing 
neutral fractions: *hi = 0, 0.25, 0.5 (left to right), with 1 Mpc 
thick slices through the reionization fields overlaid over the ob¬ 
served LAEs (green corresponds to neutral patches). Rows cor¬ 
respond to the different L“ tr -O Mh prescriptions and different 
reionization morphologies, as indicated in the labels. As described 
above, all mock surveys are shown at a fixed observed number den¬ 
sity of LAE, wlae = 4.1 x 10 -4 Mpc -4 (Ouchi et al. 2010]), inte¬ 
grated along the survey depth of Az ss 0.1. This requires a higher 
LAE duty cycle earlier in reionization (for a given Mh), in order to 
compensate for the additional attenuation of the LAEs during the 
EoR without changing their intrinsic clustering. 

There are several trends which can already be seen in Fig. [5] 
Firstly, there is an increase in the clustering of observed LAE in 
the early stages of reionization ([McQumn_et_a^|2007 Mesinger] 


profiles (e.g. Miralda-Escude|1998[|Mesinger & Furlanetto|2008a[ l, 
self-shielded systems have to be proximate to the galaxy in order to 


|& Furlanetto|2008b| | Jensen et al.l2013||Hutter et al.|2015| ), which 

is evident even when the maps are compared at a fixed observed 
number density. However, the clustering seems insensitive to the 
reionization morphology, at fixed xm (e.g. [McQuinn et al.|2007| l. 
This is mostly due to the fact that the narrow band surveys only 
localize the LAE within A z « 0.1 (corresponding to ~ 40 Mpc 
at z = 6.6), which is much larger than the typical scale of cosmic 
HII regions of EoR models. Therefore the information on the reion¬ 
ization morphology is lost when smoothing over the large redshift 
depth of these narrow band surveys. This can also be seen in the fig¬ 
ure by noting that the positions of the observed LAEs only weakly 
correlate with the thin, IMpc slices through the neutral fraction 
fields shown in green. 

Secondly, we see from Fig.[5]that the typical mass of halos host¬ 
ing LAEs also has a notable imprint on their clustering. Given that 
we do not understand LAEs at high-z, this poses both a challenge 
(in disentangling the clustering signature from the one imposed by 
patchy reionization), and an opportunity (in allowing us to better 


© 0000 RAS, MNRAS 000, 000-000 





















































The clustering oflAEs during reionization 7 


400 jjSMALL_Hli, x H i=0-0,: ; .=^fC 
:"M /I = 2x1O 11 W 0 K‘'-'$Vt; 

2 m'r&'j&TE Wl 

0 100 200 300 400 500 

Mpc 


0 ’ 
a 


400 

300 

200 

100 

0 . 


v.y«:b f>7 

7 SMALL H X-=0.0.; 


: • C •• • V‘ >■»’ V. . 

'pyt :■ v!".V. /•,' 


0 100 200 300 400 500 

Mpc 


0 ’ 
a 


400 

300 

200 

100 

0 


ViARGE_Hl( f H i = 0-0; : \'i^C 

:% h =2xio 11 iw 0 K : '" , '^:i^ 

• ■ ■'< • .••■ , /y.?,.' . . •••■...‘ f* " 

v • TT. ' 

s>#f*’ v' : ;7 

“• •*. ■ -,» 

H; '• .A;‘ • * '• '■»•• W.'^.Jir 

" '‘vA _•; 


0 100 200 300 400 500 

Mpc 


0 ’ 
a 


LARGE H L:00.;::v£ 

= 2x10' : ^V7"V'?7' : 
)r?T Ww 


400 v;._ 




100 

0 


A 1 ;•/■:' 

V-;; ;■ •• • v/y- 




0 100 200 300 400 500 

Mpc 



100 200 300 400 500 
Mpc 



0 100 200 300 400 500 

Mpc 

500 * y * jyy ;&• i vsr 

400 $£hi - 1‘ .H>jL 



0 100 200 300 400 500 

Mpc 



Jjg 

r>.. ... 


0 100 200 300 400 500 

Mpc 



Figure 5. Distribution of the observed LAEs in our mock surveys at different ionization fractions (Shi = 0, Shi = 0.25, Shi = 0.50 from left to right); 
the panels show different reionization models and L™ tr M E prescriptions: SMALL_HII, M = 2 X 10 11 AEq; SMALL_HII, M = 2 X 10 10 Mq; 
LAKGEJIII. M = 2 X 1 0 11 Afg; LARGE.HII, M = 2 X 1 0 1 0 Mq (top to bottom). All of the surveys have the same observed number density of LAEs 
n LAE = 4.1 X lO" 4 Mpc' 3 . 
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8 Sobacchi & Mesinger 


understand Lya emission from galaxies using the large number 
statistics from upcoming surveys). 

To sum up, the observed clustering at a fixed LAE number den¬ 
sity should be strongly impacted by: (i) the typical halo masses 
hosting LAEs, Mh; and (ii) the neutral fraction of the IGM, xm- 
The clustering is far less affected by the reionization model (i.e. 
morphology of reionization). In the next section we quantify these 
statements. 

3.2 Angular correlation functions 

In Fig. [6] we show the angular correlation functions corresponding 
to our models. The x-axis show the comoving projected separation; 
for reference 10 Mpc corresponds to an angular separation of ~ 3 
arcmin. The impact of the different reionization morphologies is 
illustrated by comparing the left and right columns, while in the 
rows we show the three intrinsic L™ 1 T4 Mh relations. As al¬ 
ready noted, the extreme “Most massive halos” model with a duty 
cycle of unity is inconsistent with current observations even post¬ 
reionization. The most massive halos cluster too strongly compared 
with the observed LAEs. 

The other two intrinsic clustering models (bottom two rows) are 
consistent with the |Ouchi et al.| ( [2010[ t observations, for a mostly 
ionized Universe. The correlation functions are boosted by a factor 
of ~ three when the Universe is half ionized (compared with post¬ 
reionization). On the other hand, the morphology of reionization 
at fixed xm (c.f. left and right panels) has a small impact on the 
clustering of LAEs, with the correlation function changing at the 
level of ~ 10%. 

It is interesting to notice that a higher iehi shifts the correla¬ 
tion function without changing its shapj^°1(e.g. Jensen et al. 2013). 
Thus the correlation function is robustly described by its value at a 
given scale; below we make use of this and adopt for our observed 
statistic the value of the ACF at 10 Mpc, wio. 


not change, despite a significant decrease in the l-cr uncertainty. 
For simplicity and clarity of presentation, we do not include the 
sample variance of the reionization morphology in observational 
uncertainties, which depends on the given EoR model ( |Taylor &| 
|Lidz|2014f. This EoR sample variance should be negligible in the 
upcoming wide-held HSC survey ( |Jensen et al.|2014| l. 

Figure [7] can be interpreted as follows. The blue shaded region 
corresponds to the theoretically-allowed range of ACF values; the 
horizontal orange strip corresponds to the observed value. The re¬ 
gion of overlap represents the currently available constraints from 
LAE clustering. We can immediately see that even the current ob¬ 
servations are very constraining on both reionization and the aver¬ 
age host halos of LAEs. 

Even using our extremely conservative intrinsic clustering 
model (“Least massive halos” with Mh ~ 3 x 10 9 Mq) we see 
that current 1 -o constraints from Subaru imply a mostly ionized 
Universe at z = 6.6: xm < 0.5 (see also| McQuinn et ah||2007| 
and |Ouchi et al.|2010[ l. Assuming that the mean value of the ACF 
remains unchanged, this limit should tighten to xm <0.3 with the 
upcoming HSC Ultra Deep Survey. Note that the less-conservative, 
intermediate model for the intrinsic clustering, Mh « 2 x 1O 1O M0 
already implies xhi < 0.2. 

Clustering measurements also strongly constrain the typical host 
halos of LAEs. Even assuming a fully-ionized Universe, the host 
halo masses of LAEs must be Mh < fewxl0 10 Mg, in order to 
be consistent with the Subaru measurements at lcr. The resulting 
LAE duty cycles, obtained by matching the observed number den¬ 
sities, must be at the percent level or less. These values can be com¬ 
pared with the analogous ones for z ~ 7 LBGs, which have Mh ~ 
fewxlO 11 M 0 and duty cycles close to unity (as computed from 
their two-point correlation functions; Barone-Nugent et al .|2014) . 
The disparity between these values motivates the intriguing pos¬ 
sibility that z « 7 LAEs and LBGs are dominated by different 
classes of galaxies. 


Jensen et al. 


2013 


3.3 Joint constraints on reionization and host halo properties 

In Figure [7] we show the main result of this work: the relation be¬ 
tween the average neutral fraction aim and the angular correlation 
function at R = 10 Mpc, wio. Solid (dashed) blue curves corre¬ 
spond to the SMALLJTII (LARGE_HII) EoR models. The hori¬ 
zontal orange strip corresponds to the 1 -o uncertainty from Ouchi 
et al. (2010; L™ rl = 2.5 x 10 42 erg s -1 )[^1 The width of the 
solid red horizontal line corresponds to the expected l-cr uncer¬ 
tainty from the planned Ultra Deep campaign with the Subaru Hy¬ 
per Suprime-Cam (HSC), which will probe a comoving volume of 
~ 3 x 10 6 Mpc 3 at z rts 6.6 at the same sensitivity (M. Ouchi, 
private communication). The exact position of this line is still to be 
determined: in the figure, to highlight the improvement in the deter¬ 
mination of Xm from clustering measurements by HSC, we have 
(somewhat arbitrarily) assumed that the mean value of ujio does 


10 The real-space correlation function could show a change in slope if the 
bubble size distribution is sharply peaked (e.g. |Furlanetto et al.[2006) . How¬ 
ever in practice the width of the cosmic HII region PDF combined with the 
relatively wide redshift bins, Az rs 0.1 of narrow-band LAE surveys smear 
out such a feature. 

11 These observations assume that the LAE sample is not contaminated by 
lower redshift interlopers, as argued by |Ouchi et al.|(2010f ; however, if the 
fraction of contamina nts i s f c , the ACF is artificially suppressed by a factor 
(1 — f c ) 2 - In Figure [A4| we show that a reasonable choice for f c = 0.15 
does not affect our conclusions. 


3.3.1 Implications for the relation between LAEs and LBGs at 
a « 7 

Because the HMFs and LAE LFs are steep, the bulk of the cluster¬ 
ing signal will be sourced by galaxies close to the detection thresh¬ 
olds. Therefore the simplest explanation of the different cluster¬ 
ing properties of LAEs and LBGs, is that the LAEs with « 
3 x 10 42 erg s^ 1 correspond to UV magnitudes significantly fainter 
than current Hubble limits of Muv ~ —18. The ratio of Lya and 
UV luminosities is related via the rest frame EW. Using a proba¬ 
bility distribution of EW given a UV magnitude, P{EW\M\jy), 
empirically-derived from a large sample of z ~ 3 LBGs, |Gronke| 


et al. 


(20151 estimate that 


3x 10 42 ergs -1 


galaxies should 

on average have UV magnitudes of Muv ~ —18 at 2 ~ 6, mak¬ 
ing them detectable with current Hubble surveys (see their Fig. 3; 
see also |Dijkstra & Wyithe|2012) . These empirical trends are con¬ 
sistent with theoretically motivated relations between a constant 


nicutt et al. 

1998 

Madau et al. 1998, Osterbrock & Ferland 2006; 

Garel et al. 

2015 

i. Thus one could expect that the mass estimates 


of LAEs and LBGs host halos should be more or less comparable, 
driven by the clustering properties of galaxies not much brighter 
than L™ ln « 3 x 10 42 erg s _1 and Muv = —18. Instead, the 
observed z tts 7 LBGs have over an order of magnitude higher 
average host halo masses and duty cycles than the Subaru LAEs 
( jOuchi et al.|20 1 ()[|Barone-Nugent et al.|2014| . 

One explanation to this apparent discrepancy is if the narrow- 
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Figure 6. Angular correlation functions of the LAEs in out mock survey at z = 6.6 and iilae = 4.1 X 10 -4 Mpc -3 , assuming Shi = 0, Shi = 0.25 and 
Shi = 0.5 {bottom to top in each panel). The impact of the different reionization morphologies is illustrated by comparing the left and right columns, while 
in the rows we show the three intrinsic L l ™ tT -o relations. Points correspond to measurements from Ouchi et al. (2010), with L™ ln = 2.5 X 10 42 erg 


band selected LAEs are dominated by a different class of objects 
than those selected in broad-band dropout surveys. Although EW 
distributions have only been constructed from a small, bright sub¬ 
sample of Subaru LAEs i jKashikawa et al.|2011} , they indeed ap¬ 
pear to be much flatter than those constructed from color-selected 
galaxies (e.g. |Dijkstra & Wyithe|2012[|Schenker et al.|2QI4j . 

A physical motivation for these trends could be that the narrow- 
band selected LAEs are dominated by a population of young, star- 
burst galaxies, residing in less massive halos. Since the Lya lumi¬ 
nosity of a galaxy is more sensitive to a younger stellar population 


than its 1500 A UV luminosity, bursty-galaxies have higher EWs. 
Indeed, a bursty SFR seems needed to explain the observed popu¬ 
lation of high EW galaxies at 2 ~ 3 in the models of |Garel et al.| 
( |2015] l. Spectroscopic and photometric UV follow-up of the current 
and upcoming Subaru fields could help shed light on these trends. 


4 CONCLUSIONS 

Recent observations suggest that both the fraction of drop-out 
galaxies with strong Lya emission and narrow-band selected LAEs 
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Figure 7. Relation between the average neutral fraction xhi and the angular correlation function at 10 Mpc, turn, with an observed LAE number density of 
n LAE = 4.1 X 10 -4 Mpc -3 and L™ ln = 2.5 X 10 42 erg s _1 at z ~ 6.6. Solid (dashed) blue curves correspond to the SMALL_HII (LARGE_HII) EoR 
models. As stated in the text, the IGM absorption is evaluated At) = 200 km s —1 redward of the systemic redshift. The horizontal orange strip corresponds 
to the l-<j uncertainty from Ouchi et al. (2010; L™ ln = 2.5 X 10 42 erg s~ 1 ). The width of the solid red horizontal line corresponds to the expected 1- 
a uncertainty from the planned Ultra Deep campaign with the HSC. This figure can be interpreted as follows. The blue shaded region corresponds to the 
theoretically-allowed range of ACF values; the horizontal strip/line corresponds to the observed/estimated value from Suprime-Cam/Hyper Suprime-Cam. 
The region of overlap represents the available constraints from LAE clustering. 


show a drop in these populations at 2 >6, reversing empirical 
trends from lower redshifts. These trends are qualitatively consis¬ 
tent with expectations from the EoR. However, quantifying the cor¬ 
responding EoR constrains is difficult due to our ignorance of the 
intrinsic Lya emission, and small-number statistics in galaxy sam¬ 
ples (e.g. Mesinger et al. 2015| and references therein). 

Motivated by the upcoming wide-field survey with the HSC, 
here we focus on a complementary EoR probe: the clustering of 
LAEs. The most important advantage in using the clustering of 
LAE is that the intrinsic ACFs of dark matter halos are well un¬ 
derstood, and only vary by factors of few over the range of possi¬ 
ble host halo masses. This simplifies disentangling the uncertainties 
from the intrinsic emission and the EoR, allowing us to draw robust 
conclusions from a single redshift LAE survey. 

We confirm (e.g. Jensen et a l.}2014fr that the ACF at a fixed ob¬ 
served LAE number density and xm is extremely insensitive to 
the EoR morphology; distinguishing between different EoR mod¬ 
els would therefore require more accurate redshift determinations 
with spectroscopy. We illustrate this possibility with Fig. |A5| in the 
Appendix. 


Exploring a wide range of possible host halo masses, we con¬ 
clude that current measurements of the 2 = 6.6 ACFs (jOuchi et al.) 


2010 1 imply xhi <0.5 (1-cr). Although in qualitative tension with 
the Lya fraction measurements which favor a neutral Universe, 


these two complementary probes are consistent at 1 -a (e.g. Caru-| 


|ana et al,]|2014| |Dijkstra et al.||2014| |Mesinger et al.|[20l5 i. The 

upcoming Ultra Deep campaign with the HCS could improve on 
these constraints by tens of percent, constraining xhi < 0.3 (1-ct) 
if the mean value of the ACF remains unchanged. 


Since the value of the LAE ACF is relatively low, it is possible 
to put an upper limit on the masses of their host halos (the EoR can 
only increase the value of the observed ACFs). We find that M h < 
10 1o Mq, an order of magnitude lower than the analogous value 
for color-selected, LBGs at 2 s=s 7. Combined with their observed 
number densities, this implies a very low duty cycle < few per cent. 
We suggest that this discrepancy could be due to the LAEs selection 
technique, preferentially selecting a population of young, star-burst 
galaxies, residing in less massive halos. 


We thank M. Dijkstra and M. Trenti for useful comments on an 
earlier version of this manuscript. We are extremely grateful to M. 
Ouchi for providing us estimates for the upcoming HCS survey. 
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APPENDIX A: SENSITIVITY OF THE ANGULAR 
CORRELATION FUNCTIONS TO MODEL PARAMETERS 


In this appendix, we show some alternate versions of Fig. [7] varying 
other fiducial parameter choices. In Fig. |Al| we vary the velocity 
offs et fr om the systemic redshift assuming Av = 100 km s _ , in 
we assume a L\T' T oc M 2 ' 7 ' 3 scaling for the “Least massive 


A2 


Fig-_ 

halos” curve, and in Fig. A3 we use Mh = 6 x 10 S M 0 for the 


“Least massive halos” curve. The similarity of these figures and 
Fig. □ demonstrates that our conclusions are not notably affected 
by these parameters. 

In Fig. |A4| we instead change the horizontal strip in order to 
show the impact of / c = 15% low-redshift interlopers in the Subaru 
field. The maximum impact of interlopers (assuming a random dis¬ 
tribution) would be an artificial su ppressi on of the ACF by a factor 
(1 — /c) 2 - Although Ouchi et al. (2010) argue for f c ~ 0%, and 


find no evidence of contaminants in ongoing spectroscopic follow¬ 
up, they quote a strict upper limit of f c < 30%. In the case of in¬ 
terloper contamination, the allowed range of Mh values broadens 
somewhat, but our main conclusions remain unchanged. 

Finally, in Fig. |A5| we show how LAEs can probe reionization 
morphology, provided they are better localized in redshift space. In 
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Figure A2. Same as Fig. but assuming a cx ^ 3 scaling for 

the “Least massive halos” curve. The similarity of this figure and Fig. |7] 
demonstrates that our results at fixed M^ are not affected by the precise 
mass scaling of the Mh relation, since the clustering signature is 

dominated by halos with masses close to Mh- We do note that these halo 
mass-to-luminosity scalings, if extrapolated to the bright end of the LFs, 
would result in very different predictions (see Fig.QJ. 
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Figure A3. Same as Fig. [7] but Mh = 6 x 1O 8 M0 for the “Least mas¬ 
sive halos” curve. The similarity of this figure and Fig. |7] demonstrates 
that our results are marginally affected by a further decrease of Mh below 
3 x 1O 9 M 0 . 


X HI 

Figure A4. Same as Fig. [7] but assuming a fraction f c = 0.15 of randomly- 
distributed, lower redshift interlopers when plotting the horizontal strip cor¬ 
responding to Subaru observations. Although the allowed range of Mh val- 
ues broadens somewhat, the main conclusions remain unchanged. 
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Figure A5. Predictions for cjio as in Fig. [7] but assuming a survey depth 
of A.z = 0.014 (instead of A z = 0.1), as could be achieved with Lya 
spectroscopy. The dependence of tojo on the EoR morphology is greater. 


this figure, we take a survey depth of 5 Mpc (A z = 0.014 instead 
of A 2 = 0.1). The curves for the two EoR models diverge much 
more than in the previous figures. This survey depth (correspond¬ 
ing to An ~ 600 km s _1 in the rest frame) is broader than the ob¬ 
served spread in the systemic line offsets (e.g.|Shibuya et al.|2014| 
|Stark et ah 2014, Sobral et al. 2015); hence such measurements 
can be achieved without metal line detections. These results illus¬ 
trate that wide-field LAE spectroscopy, either through follow-up 
or through dedicated instruments such as the Subaru Prime Focus 
Spectrograph, could be very powerful in discriminating between 
EoR models. 


© 0000 RAS, MNRAS 000, 000-000 



















